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ABSTRACT 

Fast voltage transients have been found capable of disrupting the 
operation of a proposed solid-state, ac-dce inverter system. A method is 
proposed for studying the propagation of such transients in an involved 
d-c power distribution system connected to such an inverter. The use of 
computer simulation of the system is proposed and detailed. Mathemati- 
cal equivalents of various circuit elements are developed, and signal 
flow graphs are used to write the system's equations in a form adapted 
to the computer. Use of this method is shown to allow a large number of 
system locations to be chosen as inputs or as points of measurement, 


even when the simulation is programmed for a digital computer. 
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be Introduction. 

1.1 The Bureau of Ships, alert to the potential of static power conver- 
sion equipments for use between d-c and a-c power distribution systems, 
has conducted studies and has contracted for studies of the feasibility 

of shipboard installation of such equipments. The advantages of static 
power conversion devices over motor-generators are numerous. "The static 
inverter, being solid state throughout, has no moving parts, needs no lubri- 
cation, is lighter, uses less floor space and has greater efficiency than 
its rotating counterpart." [5] The studies initiated by the Bureau of Ships 
have included the factors of cost, reliability, maintainability, space, 
weight, and environmental considerations, as well as a-c and d-c systems 


performance. 


1.2 The purpose of this paper is to plan in detail the study of one of 
the problems that has been found in pilot systems where such equipments 
have been installed. That problem is the presence of fast, high voltage 
transients at the d-c input terminals of the inverter caused by faults 


and fault clearances. 


1.3 Although the study described herein is concerned primarily with the 
d-c distribution system and only remotely with the conversion equipment, 

a brief description of the sensitive solid-state device is necessary to 

the understanding of why the voltage transients are of concern. The 
principle load-carrying element of the solid-state, static inverter is 

Phe silicon-controlled rectifier, SCR, a four layer device which utilizes 
three p-n junctions and three terminals to perform its functions: reverse- 
bias blocking; forward-bias blocking; and controlled switching to the con- 


duction state from the forward blocking state. High voltage transients 
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across the SCR, when they exceed the peak inverse voltage of the device, 
can appreciably shorten its life, and if sufficiently severe, can cause 
its immediate failure. E. K. Howell, of the General Electric Corporation, 
describes an additional limitation of the SCR when subjected to fast tran- 
sients, "...false triggering by the dV/dt effect. When in the blocking 
forward-biased condition, the center junction of the SCR is essentially 
a capacitor connected between the anode and the gate. If the anode volt- 
age is changed, a displacement current will flow through this capacitor 
equal in value to the product of the capacitance and the rate of change 
of anode voltage, dV/dt. A high dV/dt caused by a fast transient can 
trigger the SCR into conduction" [2] - Such spurious gating of an SCR 
at the d-c side of an inverter may result in one or more of the follow- 
ing effects at the aec output: 

Frequency drift 

Voltage variation 


Harmonic distortion of the a-c voltage waveform. 


1.4 Specifically, the problem was to determine the magnitude and rate of 
rise of voltage produced at the inverter by sudden interruption or initia- 
tion of large currents elsewhere in the d-c system and to determine what 
the effects would be if the system were slightly modified. System simula- 
tion with a computer was the method of analysis chosen. The bulk of the 
study is adaptable to either a digital or analog computer. Only in the 
final stages, section 5.1.2., is a commitment made in order to solve the 
particular problems of this particular system. The advantages and dis- 
advantages of the chosen course, digital computer simulation, are discuss- 


ed and details of a suitable program (Fortran '60) are outlined. 








ne D-C Distribution System. 


2.1 General Description 


A block diagram of the d-c system showing major components and inter- 
connecting cable is presented in Fig. 2-1. The basis for this diagram are 
references [3] and [4] - Not all loads indicated in those references are 
included in Fig. 1, nor are they considered in any of the subsequent develop- 
ment. Their omission is justified by the following considerations: 

a. The emphasis of this investigation is upon developing tech- 
niques and procedures for studying fast transients. Given the distribution 
system of a particular class of submarines, and using the suggested cir- 
cuit models and simulation procedures, the response of the system to a steep 
wavefront current can be readily predicted. The validity of the techniques 
is not dependent upon the number of loads chosen for this development, but 
rather upon the accuracy of the component circuit models, and upon the 
assumption that all significant dec loads and transmission components can 
be represented by the proposed models; 

b. Omitted loads are in such a location that they would not 
affect the transients at the chosen points of interest were they included; 

c. Circuit protection equipment (fuses, circuit breakers, low 
voltage relays, etc.) are omitted. The transients of concern are very 
fast compared to the reaction time of these devices; that is, the trangs- 
ient will have reached peak magnitude before the protective equipment can 
operate to interrupt the transmission in any branch. 

d. The detailed plans for those submarines for which this study 
may be of interest are classified. 

From the block diagram it can be seen that there are two distribution 


system, port and starboard, tied together electrically at the battery. 











laure, 
eee WS ZI-VISS 


*uleysAs uoTyNgTW4STq o-p 
au} Jo weisetq yooTg 


T-¢ eingty 
(hOS9-V9OSS 


(DO0O08-V9SS (Q)059-VOSSs 





WOd-—-Sds (ZI008-V4OSS I(€)009-VDSS 







SOLS 
uoLs(ndody 
iC 


£49129 - 


P9246 -SdsS 1(22009-V49SS [(HO08-V9SS 





() OSS -V¥ IS 






JoyeuaUag 


Jasaiqg'étow 2 Jaury 


CDS ley 26 Ge ESIC 








With battery breakers open, each system can be operated independently 
with power supplied by the a-c to d-c conversion equipment (see section 
22h). 

As is standard for all shipboard electrical installations, the d-c 
system floats with respect to the ship's structure, and thus comprises a 


complete electric circuit. 


2.2 System Components 


2.2.1 Battery 


The battery provides a high capacity energy source with nominal out- 
put voltage of 250 volts. Its primary functions are: 

a. Provide power for all control and instrumentation circuits 
during reactor startup; 

b. Provide power to essential circuits, a-c and d-c, when 
reactor has been "scrammed" or is otherwise shut down, or is functioning 
at too low a power level to operate these circuits. 

The linear circuit model of the battery is developed in Appendix C. 
Note: Considerable use is made of the Laplace transform wherever the 
occurrence of linear relationships among circuit variables render its use 
advantageous. Its use, however, is nowhere essential; it can be replaced 
completely with analysis in the time domain. The use of upper-case letters 


for system variables indicates the transform of those variables. 


2.2.2 Trim and Drain Pumps. 


These pumps, one fed by the port bus, the other by the starboard bus, 
are representative of any of the d-c motor loads. They are the largest 
motors, in terms of both horse power and mechanical inertia, that are 


frequently on the line. Their linear circuit model is also developed in 








Appendix C. 


2.2,3 Emergency Propulsion Motor Feeder. 


This feeder is included to show the effects of any load feeder which 
is open-circuited at the load end. The circuit model for this feeder is 


identical to that of all the other transmission cable discussed in 2.2.5. 


2.2.4 Static Power Suppl SPS). 

On nuclear powered submarines currently in commission, power conver- 
sion between the a-c and d-c systems is effected through motor-generators. 
There are usually two such machines, each rated at 300KW, and each cap- 
able of two-way operation. In this study each motor-generator is consid- 
ered replaced by its static equivalent, the SPS. It is the solid state 
devices which comprise the load-carrying circuits of the SPS that are 
sensitive to the fast transients. The passive circuit model of the SPS 
depends upon which of several possible circuits is used. In Appendix C, 
one such circuit has been shown, and the branch termination representing 
the input impedance of theSPS has been established from this circuit, 


based upon the selected mode of operation. 


2.2.5 D-C Bus and Load Feeders. 

The d-c buses, port and starboard, consist of single-conductor, shield- 
ed, copper cables (SSGA), the number and size of which in any part of the 
system depend upon the rated load of that part. Between the distribu- 
tion panels and individual loads, two-conductor, shielded cables (DSGA) 
are used. In Fig. 2-1 the cable type designation is followed by the 
nominal size of the cable (expressed in thousands of circular mils); the 


number of cable pairs is shown in parenthesis following the size. 








The transients of concern are short duration, high frequency voltage 
components superimposed upon the normal d-c line voltage, or upon slower 
transients. This means that cables must be considered as transmission 
lines. In Appendix B, expressions for the distributed parameters of the 
cables are developed, and in Section 3, a lumped parameter equivalent 


model is established. 





3. Modeling the Connecting Cables as Transmission Lines 


3.1 Transmission Line Simulation. 

The sets of cables were reduced to a single transmission line by re- 
placing the distributed parameters of the individual cables with equiva- 
lent distributed values, and by saying that all losses and interactions 
occur only in one conductor of the pair. On this conductor, all voltages 
were measured above a ground which is designated as the electrical mid- 
point of the real, balanced circuit. 

The duration of the transients considered was in the range of one to 
fifty microseconds [6] - Since this is the same order of magnitude as the 
time required for a wave front to propogate the length of the cables, the 
propogation time must be taken into account in the model and the cables 
must be considered as transmission lines. 

A transfer function representing these lines is derived in section 
3.2. This function is shown to be directly frequency sensitive, and to 
contain some parameters which also are functions of frequency. The fre- 
quency sensitivity of the transmission lines and their representative 
lumped networks must be emphasized. A transient can be represented as a 
wave form or as a frequency specturm. A transient passing down a trans- 
mission line has its wave shape modified; the change can be analysed in 
terms of the transient's spectrum and the line's frequency response. Two 
devices with similar frequency response characteristics will have similar 
effects on transients passing through them. The purpose of this model is 
to duplicate the transient response of the large distribution system, so 
its components must have frequency responses matching those of the large 
system's components. 

In section 3.3 an equivalent circuit is developed to duplicate the 
line's transfer function. With three elements, however, the equivalent 


8 








circuit is correct for only one frequency. One hundred kilocycles was the 
frequency chosen since all the lines could be simulated at this frequency 
and be approximate simulations at lower frequencies. If the results, based 
on the 100 kc model, indicate that higher frequencies may be important, 
the range of the model can be slightly extended, to about 500 kc. If 
still higher frequencies are required, the form of the equivalent circuit 
must be changed, and it will either be valid for high frequencies only, 
or will require extra elements. The effects of picking the 100 kc model 
are illustrated in section 3.4. 

The equivalent circuits replace the cables in the system model shown 


in section 2.1. 


1 


3.2 Transmission Line Equations and Transfer Functions, 





3.2.1. When a voltage, Eog(t) » is applied to the input terminals of a 
transmission line and a current is sent into the line, the voltage and re- 


sultant current are initially related by: 
E <x) = Lix)+ZLo (3-1) 


This is one of the solutions of "The telegraphers’ equations": 

dE /ox =-(L-R+L-d1/dt) 
OT/dx = -(E-G+C-3E/dt) 

Here R, L, C, and G are the distributed parameters of the transmission 
line, the values seen looking at a very short section of the Line and 
considering it in relation to the relatively long, straight lines joined 
to it. Their units, typically, are: 

R: ohms/meter L: henrys/meter 


G: mhos/meter C: farads/meter 


l 
This section is based on material from | 1.2] and 
9 








Taking X as the length down the line from the input terminals, 
—ObX 
E(x) = E@)-€ = El) exp(-Tx) (3-2) 
ae and § are the transmission line constants, whose values depend on 
the lines’ geometry and materials. ~ Wa is the line's surge impedance or 
characteristic impedance; 6 is its propogation constant. The exponen- 
tial term, in 3-2, is negative because the direction of power flow is in 


the direction of positive X  , away from the source. 


Decale If the applied voltage is taken to be sinusoidal, of the form 


Eg x p(jwt) and if G is negligible, Zz and & are given by: 


a (EEA 423 
Ms % ne ohms (3-3) 


¥ =\(R+joL)jwCc) 








(3-4) 
ae as =(jwC)Z,) 

a lL Le yee 

= al 7 ) (FS + wCZ,)| 

= a+jB 

a _R 
Q= ZTILp (3-5-a) 
a ole, + a 


If Go is large (it is 30 or greater in the system being modeled)! , 


2 
Appendix 8.4 shows their calculation. 


1 
Table B.3, Appendix B. 
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Zox = (3-5-c) 
B=wYLC (3-5-4) 


Now equation 3-2 can be written: 


E(x,t) =|EotexpGwt-jwic x -ax)  G-6a) 
and, 


|E(o,t)| = 


(3-6-b) 








ve Zz, Ste g(t) 
This shows that at time t, the voltage at X lags the voltage at X 
equal zero with phase angle equal to (IK , and the magnitude has 
been attenuated in traveling from 0 to X_ by Clee ly.) The same 
happens to ALIGGE 

If the input is other thansinusoidal, or if transients are of in- 
terest and Laplace transforms are going to be used: 

Sf 
5 a = ee 
With 
Zhe, t= arsliC 

Now, in transformed form, equation 3-2 is: 

E(x,s)= E@s) expC-ax)-exp(-sx\LC) 

Again, expCax) represents attenuation. The second exponential 
term now represents a time delay, showing there is a time lag equal to 
yLc before an event at X = 0 has an effect at distance xX . 

A new quantity, Vo, the phase velocity is easily seen in this intere- 


pretation: 


eee | 
Ve yc (3-6-d) 
1 


1 








This is the speed with which a wave travels down the line, and: 





Sx 
E (SS) =-E@s) exp Gare i=) (3-7) 
In these, E (0S) is the transform of —<o— Eg (t) 
Zona 
Breede Up to this point, a single wave of voltage and current has been 


pictured, traveling in one direction. The conditions governing it neuen 
characteristic impedance and the propogation constant. Let this be called 
the forward wave, and subscript it with an f. At the end of the line 
there is, in general, a loddeuien some impedance, ae At this load, 

to the right of the f-cerminals in Figure 3-1, this impedance defines the 
relation between voltage and current. But to the left of these terminals, 


Za6 relates them. 


si Lt) f 
Z9° e ¥ + = 


' Ty 
al “Boat ah) 
| 
SOURCE , LINE ” LOAD 
— | 
Figure 3-1 


Transmission Line In a Circuit 
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Loosely speaking, if Ly is larger than Jone Ty Zo is larger 
than Ep ,» and an increase in voltage begins to propogate as a wave 
(with speed VY ) from the load toward the generator. Simultaneously, 
Eep/Zy is smaller than HES. and a decrease in current is required, 
which propogates as a negative current wave. Conversely, if ya is less 
than Lowe Ip is too small for equilibrium and an increase in current 
propogates as a wave of Re ieive current accompanied by a wave of negative 
voltage. 
Speaking more precisely, at any time: 
E,=1yZ, = Eo, : 
TE ties ES (3-8) 
Er(Q)=Le(0):Z, 

These requirements are met by adding the returning, or reflected, 


voltage and current waves, Eo), T,), 0 the forward ones. The ree 


flected waves are found from the load's coefficient of reflection, ky : 
kp = Z¥-ZLo (magnitude =) KS RY Zt 0.) 
Z4+ZLo 


Er(9) = ko Ee (8) 
Lr) = —ky Lp) 35) 
Ge E p(x) +E OX) 


Tm= Leary) 


The quantities determined by equations 3.9 satisfy the requirements 
of 3-8. The return waves propogate toward the generator with the same 


attenuation and delay per unit length but with: 


Er(xyy =-4Z: 1, &%) (3-10) 
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The distance traveled from to X is J. » 80 (q-x) replaces X 
in the exponent. 


Taking the sinusoidal case, 


Fr (xt) = Ey(2,t)-exp(-1(2-») 


ke ° HAG ks KO(-¥CL- x) 


Eo. 1p a 
kg: bya 7. S, ; exp|j (wt-B2d- -x)} ‘expFa(2k- x) 


showing the phase difference ie Rees the applied signal and the returning 





wave at any time, e » and distance, X. Similarly, 
E-(X5) = ky -E-(0,5)xp|-a(2 2-x)- cn (3-11) 


Note that if the applied voltage, E4(t) | had included a switch or 
a unit step function, LU » the time expression found by inverting 
3-11 would include a unit step, U(t-24=%), 


With time, t, less than 20/r , at X on the line there is: 


Ex) =Evosyfexp(-xCa+&)) + ky exp(-@et-vlas wl 
iG) = EON expCacas S) £))- ky: expl- eta &)] 


The reflected wave reaches the input, ie terminals at (ee 2a. and 
again there is, in general, a mismatch and a coefficient of reflection, 

ke » must be defined for the generator terminals. (A match occurs 

when the terminal impedance, Lo or Ly is equal to Flags and the reflec- 


tion coefficient, ke or ky is zero.) 


Ka = Zaq- Lo (3-13) 


ZotZo 


14 


Now there is an additional wave traveling from the generator towards the 


load, which is subscripted ft while the original is fh 


L ¢, (0,5) = kg Eales) = ka Ky E ¢, (0,5): exp G2k(a+ +) 


(3-14) 


Tp, (0,3)= “kg Ty, (08) = kg Ky Tp, (0,s)-exp-2 hla +3) 


| SQ ae eos 
With t less than ar and = At aF; 


Li) = Ep (XS) + Ep (x5) + Ep, (%5) (3-15) 


=E (05) @x p(-X) +e exper t-x)) + key ky ex pfs (2 Q+%) 


ILES)= EOS) expcry)-ly exp Ca (ak-x)+kg kyexp-r(20+%) 


At the L terminals, the current entering the line is the algebraic sum of 
the current waves and the terminal voltage is the algebraic sum of the 
voltage waves. To the left of these terminals, all conventional lumped, 
circuit conditions are fulfilled. 

Before the notation becomes overwhelming, the conventions in use 
thouls be reviewed. An ideal voltage source (zero internal impedance), 
Eg(t), is connected through a generator impedance, La , to the L ter=- 
minals. (Ce) , equal to 722Falt) is the fraction of that voltage 
seen at the terminals; E94 (09) i E.(0$) are the transforms of these 
quantities. E (03) is the first forward wave, E 5&5) at X= 0. E;(s) 
is the sum of all the voltage waves at X= 0 and JT; (s) is the sum of 
the current waves there. Likewise, Ep(s) and Ty(s) are the sums at 
x=X, 


15 








Any change in Eq affects Ey only through its effect on Ef, until 
the changed Ep, wave has gone to the end of the line and returned as eae 

The fe wave, reflected at the generator by the V| wave, eventually 
generates an V2 wave, and the process goes on forever, with ever increas- 
ing attenuation. 

At a time when there are an even number of waves present on the line 
(every forward wave Can be matched to a returning one), the line voltage 
is: 


Exs) =E(o sfexplr) ky exp(-#(2.2-x)) + q@kgexplo(2h+)). . 


(3- -16) 
rp hy expea(ankerye ky” KP expCrlenkezt-2) | 
(2nsi)*h Term (2n+2)th Term 


The terms can be grouped together as: 


ED = ELCK Oe aS) — }Codd terms) + + )(even terms) 


(3-17a) 


with t ¢(x,$) = E(os)| expCsa)+ ka kgexpCralin))s. : 
* kis ka exp(-r(ankex) | (3-17b) 


Each of the terms in 3-17b includes a delay which, when the equation is 


inverted into the time domain, makes that term equal to zero until the 


andax 
We 


time is equal to or greater than Therefore "1 may be allow- 


ed to grow very large, and as it does, 


ex p(-¥%) tice 


LE 6(xs) a | - Kg kg Exp (-24 2) ii 
Similarly, 
— _ kpexpCr(2a-2) 
Er (XS) = Tea dey Exp C2Y 2) a 
Also: Lotx, 3) = os S) Tp(x5)= 5 (3-19) 
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From these, 


oe expr) ky exp-s(22-x) . 


I-kg kg exp(-2yX) E (0,5) (3-20) 
_ _@xpCrx)-kpexpCsc2k-X)  E¢o,8) 
10s) = I~ ky kg Exp (242) ae 


3.2.4 Setting X = 0 or to Q gives the values at the input or load: 


| § Ba kg exp(-27 2) 
E:¢3) = E¢o,3) = Ka ky OXp(-2¥Q) 


(3-21) 


Ecos), \-ky-expC2rQ) 
Zo —I- kgke exp(-2 0) 


7 ex p(-72) (1- kg) 
i) ill COD Sep oy 


- Eros), expGs&)(\- Ky) | 
146s) Zo oy ka Kp @xp(-2¥ 2) 


Lics) 


Rearranging these into a more convenient form, a voltage transfer function, 


a current transfer function, and an input impedance are defined: 
Ey(s)_ | +k ~ 22) 
Ay = Eis) exp(rh) b K Exp x0) ~ Ly fas WZ sinh vf 
=. 1G) 
= a= = ExNGY)- nto ) cosh +ZgSinh EX (3-23) 


Z, = ErG)- Zy+ Ztanhrl 
; Eg(s) Zo+ Zytanhvs (3-24) 
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Any two of these equations define the system to the right of the terminals 


in Fig. 3-1. 


3.2.5 The trigonometric functions of jC or j¢ are now substituted for 
the hyperbolic functions of Y , using the following approximations, 
which are discussed in section 3.4.2. 


With small & :; 
Sinh@+j@AlL = Jj sin@l = jf sin(Y%r) 


teni(x+jP)%e = j tan(P4/) = jtan(*&er) 


The three forms of the transfer functions are now in the sinusoidal, 


steady state, form. The voltage gain, equation (3.22), is now: 


Ay =I/(COS Yr + jf SZ, Sin Yar) (3-25) 


If the load is inductive, its impedance is pul and the voltage gain 
has the form: 
Ay = (COSY + ZfwL sineYa 
This function has a pole each time tang = - wh - This occurs once in the 
0 
second and once in the third quadrant every time the electrical angle coh 
varies through Tf radians. Each time, the value of wh is closer to 
of or an (1 an odd integer). Fig. 3-2 explains this drift with in- 
, Gh. Zo 
creasing &) . The minimum magnitudes are when ta nay = Cop? and show a 
similar drift. Fig. 3-3 shows this frequency response for the case in 
£2 , . 
which “[ is larger than 1 ohm/henry. For a given f and sy and a band- 


width determined by WW, (lower) and W,,(upper), there will be m-n poles 


of A, where Mm and YH are the largest integers such that: 





ry ams ne Unk 
air aw TT 
If the load was capacitive, there would be as many poles, but in the 


other quadrants. 
18 
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Figure 3-2 


wi = 
is = NE = “Sees 
Y= St O ies = w 
| a an vt avr 3tr UE 
wo = Shy = 


Illustrating the drift of the intersections 
d ith H 
of u, an Os wit ai 





Voltage Gain With An Inductive Load 
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With a mixed load, Ly — R+ JW L and: 
R+tJjwk 
Av = Reoseb+ j(wl Bones LZ. Site GE 
As shown in Fig. 3-4, the magnitude of this function varies - 
ly between a growing maximum and a minimum which rises from zero to unity 
asymptotically. The magnitude is equal unity to every TT radians, 


and the phase angle varies clockwise from zero through —TT around to 





zero. The magnitude maxima occur when wk is slightly greater than 
~ ~1 Vota gu -| 
tan| Lars ee minima when wh is slightly less than tai) 
Figure 3-4 
Magnitude and Phase of [i With 
An R-L Load 
Pe 
eee 
—— 
— 
IAI ae 
Vv en 
| a ee eit aan, Qe = SS Se 
— ao 
Paik 
O | wh 
ims 21T 3T7T Ve 
Tei a ee ee Se ee 
PHASE 
ANGLE 
a O TT 2 ZI wl 
Vv 
Ay 
-1T a a 
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If the load is a simple ee 


(ay 


the voltage transfer function 





is: Ope A 





Fig. 3-5. 


|e 


Figure 3-5 
A, with a Resistive Load 


eee 7 = 


When the load is inductive, the current and impedance functions have 


the following form and are shown in Fig. 3-6. 


At 


\ /(cos Q2 ~ ol sip wS ) 


Lo u 


7.- j wh +Z. tan Ss 


y Zon = iD) Fan ws 
= 


These three functions have been sketched to illustrate how they vary 
with Cd). The effect of adding resistance to the load has been shown to 


dampen the function's extreme swings but not to make it less oscillatory. 


3.3 Equivalent Lumped Network 


With a steady sinusoidal input, the transfer functions are: 


y= Zy/(Zy cos BA + | Vhs stn ES 


i=Z/(Z. COO kj a) - wh) (3-26) 
a tj Zo tan a 
Zi= Lo ot | Zt aa ee 





Current 





Gain, 
Ai | 
© 
-| 
Figure 3-6 
Current Gain and Input Impedance 
With Inductive Load 
Input 
Impedance, 
Zi 





The network of Fig. 3-7 has the same transfer functions at one 


frequency, G),,if its elements are determined by the following equations 


Meco x 
Z; = jZ.tanPir Y= isipSt om 
7] 





Figure 3-7 
Equivalent Lumped Network 





When the matching frequency is chosen so that ond is less than TT, 
the series elements will be inductive and the parallel element will be 


a capacitor, with: 


coe 
ee ou), J sie re ; 
Ls 7 5 tan yar Cp ~ Oy Zo SS) 


The frequency response of such a network, connected to an inductive load, 
is shown in Fig. 3-8. When wd is in the third or fourth quadrants, 
the series elements will be capacitors and the paraller element an inductor. 
With the same load, the order of the transfer function is raised, as shown 
in Fig. 3-9. 

There is not a convenient, explicit relationship between the equivalent 
circuit's transfer functions and those of the transmission line at frequen- 
cies other than CO; . Even the poles are not usefully related. For 


instance, the input impedance pole of the L-C-L circuit (C3 in Fig. 3-8) 


can only be related to the transmission line's parameters by the equations 
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Voltage Transfer and Input Impedance of 
The L-C-L Equivalent Lumped Network 





Figure 3-9 


Voltage Transfer and Input Impedance of 
The C-L-C Equivalent Network 


“) fi. 


Le 














circuit's Ay rises monotonically. The two are equal at GW, . Over 
the lower frequency range, the two curves are bounded by two straight lines. 
It is because of these bounding values b, and b. , illustrated in 


Fig. 3-10, that the L-C-L circuit is preferred. 





C4 Try Co| W2 Ws 
ra 
Figure 3-10 
Comparison of Voltage Transfer Functions 


t = transmission line's 
€ = lumped equivalent line's 


b,-b, is line from line's dc value to W 
b2-b,is line from equivalent's dc value to line's minimum(u,4) value. 





C2= fe. tan cor 8 C4. = ae SCO oy 


| a 2 
ae 
Ws = OZ (~ cos ak) +cL Sin SE 


A more detailed analysis requires comparing numerical solutions for a 





particular line's parameters and for the equivalent circuit derived for a 
selected matching frequency. With this approach, an optimum equivalent 
could be defined and designed for a single line, but it is doubtful that 
such an optimum would be more meaningful than this procedure, picking a 
matching frequency and finding the equivalent for it, within specified 
limits. 
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3.4 Approximations. 
3.4.1 Equations (3-22) and (3-24), repeated below, give the voltage 


transfer function and input impedance based on the wave equations. 


Ay = Zl cosh ¥¥+ Z, Sinh v2) (3-22) 
7-7 2+ Z. tenn 1h 
o"Z, + Zp tanh ¥ (3-24) 
With the basic circuit of Fig. 3-1, where an ideal generator or source 
is connected through a generator impedance to the loaded transmission line, 
the load voltage, E(Q,t), can be related to the input voltage, Bea by 
the equations above and to the applied voltage, Eo(t). Alternatively, 


equation 3-21 maybe taken as the starting point and the same result de- 


rived, that (in transformed form, with Q= $f - oT; 


= “Zo _) Cit+kp) expesr ; 
EUs) ESO (ZA) (niggers) OM 


tele Ar i in sly) 1+ Kpexp C287) *.- 


+(kg ky) exp(-2 ns7T) +, , 


With the following values (those for line number 3 and with an assumed 


source and load, 
Q=67.056 m Zo= Rg=I5a Ze= Ry= 500 
L gist = 0.3634 ae Cie = 2569 za 


the values of the parameters in equation (3-28) are 


Zo =11.894n kg = 0.1159 


WV =0.3273 «105 ky = 0.615 66 


T = 2.0487% 10 © sec 
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Eg(s) — (0.7144-€xp2.0487*10%): 2 Ex p(-211+2,0487+10°S) «(0.071 ) 
Eq(s) n=O 


If the input is applied through a switch, Ea(t) = U(t), then E(t) 
is as shown in Fig. 3-lla. As M gets very large, Ey(t) approaches 
0.7687 volts. 

Using the equivalént circuit, the values are: 


L\= 11. 3B uhy C2 = 9.464xI0 °F 


and 
e 


Egy = sg 1954 O_o 
E4(S) Soe 5.711708 S+7. 6482*10'%5 + 5.3034% 108 


ee a ee 


— 3.2828x\0° a 8.8044 «10""— 3. 30x10°S 
S+4.1BBx10° —, (47-62 r108)* + (8.28% 109)" 


Again, with a switched input, 


Ep(t) = (691) + (.078)|1 - ex p(t /2.39x107)| 
—(.945)|exp (- t/i.3116* )] sin(6.28 «10% + 1.08) 


— (.399)/exp-t/i.six 10)|'sin(e.28«10" t) 


which is shown in Fig. 3-11b. 

Figure 3-11 shows the resemblance between the instantaneous response 
of the lumped element equivalent circuit and the delay-and-reflection 
response of the actual transmission line. If the equivalent circuit's 
frequency response could match the transmission line’s at higher free 
quencies, there would be even greater similarity. 
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Figure 3-11 
Comparison of Response to U(t) Input 
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The replacement of the hyperbolic functions in section 3.2.5 is 


3.4.2 
¥=A+J5B 


justified by the following. 


tanhes — Ceinh xf ) = Si +) ue 
. Ccoshsd +1) ‘2 ante 
~ ginhal-cosB+icoshad- sings 
i+ Coshad:cosBl+ j sinhak- singe 


— Vsinhtal + sint Bg /é ao tan sings ) 
coshax+cos@ 
sinh al~ al and cosh ala \ 


sinh .| 2 .l002 
cosh .-l - |.005 


For small al , 


(for example: 


so that when Sin BL Dal, G90 and ALO = jA 


In the system under consideration, the smallest value af the ratio, 


Sin Be with 100ke values, is 67.85 (Table B-3, Appendix B). 


Z 





@) 
The statement that @=90 ,» is in error by 0.84° in the worst case. 


The real part of A, equal to 1.48% of its magnitude, is being dis- 


carded. 
Similarly: 
sinh f= sinh ak. cosph+jcoshal-singl 
= inh? al + sin*el Le 
eS = tant) — tq sin Be cosh of = 90° 
o= tan io = ln OL cos@z _ 
As above: 


sink? al+ sintae lo x j singe 


As singe goes toward zero, where the approximations are poorly justified, 
the approximations go to zero also, where their exact counterparts approach 
very small, real, values. The implication is that when dec, steady state, 
conditions are considered, the equivalent-circuit model of the trans- 
mission lines would not play any part in the distribution of voltage and 
current throughout the circuit; when transient conditions are considered, 
the real lines would have a damping effect which the model doesn't show. 
The damping of the motor loads is an order of magnitude larger, however, 


and this more important effect is retained in the model. 


3.4.3. When the equivalent distributed parameters of a length of line 
have been calculated, Section 3.3 derives an equivalent lumped circuit. 
With the following manipulation, a simpler method is shown for arriving 
at the values for the lossless L-C-L circuit's elements. From Equation 
3-28: 


fe= £2 tan gis ah Cp= =~ sind 


vr 
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where Le and Ge are the series and parallel elements of the 
equivalent circuit: (CO, is the match frequency, 2 is the line's 
length, and Tae is the lines’ characteristic impedance. From 3-5-b 


and 3-6-d: Zo=\€ ) Mae ie 


where E and C are the distributed parameters. 


Therefore: 
= Before) _ aL (tanger) 


a: ( a1) (fis (3-30) 


Cn= gp (sine) =a; (sin 24) 
p = — (Bb) = © (ek) 


If the loads terminating the line had a very small resistance, so 
that the lines’ attenuation was an important part of the model, a series 


resistance could be found from equation (3-27) with tanh replacing 
Jtan ws or j tan BL ; 


At the matching frequency: 


Z,=Rstjurlg = Z, tanh 


: Pie joe = : 
leg - meme ) ¥ ={R+iaL)Ciw,o) 


SO Ve GL ia &(R+jw,L) 
? Zz. 





and Re - 9 tanh 24 ~ Rg tan we 
>a 2 2 os 2 AY 
2 2v 
If the matching frequency is taken as zero, then Las aa aL P Roe RA 


and ores 16 - A comparison of the voltage transfer function of this 
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equivalent circuit and of the transmission line is shown in Fig. 3-12. 
The pole, Wz , of the equivalent circuit has moved to the left of where 


it was in Fig. 3-10. 





| 


TRANSMISSION / 


LINE “Cy//" | EQUIVALENT 
CIRCUIT 


os 

i 
Figure 3-12 

Voltage Transfer Functions with Zero Match Frequency 








VAT 


a 





3.4.4 A basic assumption has been that the transmission line's dis- 

tributed inductance and capacitance are independent of frequency. Or, 

in deriving the delay form-equation (3-18) - it has been assumed that L 

and C were independent of ar and GE - But Appendix B shows that 
dt dt 

the distributed inductance varies with frequency because of the skin and 


proximity effects. Table B-l, summarizing this, shows that the dec in- 


ductances of the system's cables are from 1.10 to 1.14 times as large as 


the 100 kc values. The rate of decrease, - oh appears larger at low 
frequencies. 
Ze and a are the two variables which have been emphasized in this 
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4. Electrical Circuit and Mathematical Description. 


4.1 General. 

The d-c distribution system is described in Section 2; circuit models 
of the various system components are developed in Section 3 and in Ap- 
pendices B and C. The composite circuit model of the system is shown in 
Fig. 4-1. Although shown as a single line network with a common ground 

| 
return, the ground is actually the electrical midpoint of a balanced 
circuit, as noted in Section 3.1. The T-sections which are the equi- 
valents of the various cables are each designated by a single number in 
parenthesis; the number is the subscript identifying the L and C ele- 


ments that make up the T. The lower case letters on the figure are used 


for node identification. 


4.2 Mathematical Description. 
4.2.1 Requirements. 


The requirements imposed upon any model of the system, either bread- 
board or mathematical, are as follows: 

a. The model should be sufficiently flexible to allow for 
additional loads and/or variation of the parameters of both loads and 
line. 

b. Certain of the circuit nodes must be accessible, either as 
points at which the disturbance is initiated or as points of interest at 
which transients will be measured. 

c. The model should permit application of various initial 
conditions at appropriate parts of the system. 

d. If a mathematical model is to be used, it should be readily 


programmable for analog or digital computer solution. 
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4.2.2. Signal Flow Graph. 


The requirement of flexibility dictates against a breadboard model. 
A mathematical description which meets all of the above requirements and 
offers other advantages as well is phen flow graph. The addition- 
al advantages offered by the flow graph are listed below: 

a. It guards against algebraic redundancy in the formulation 
of inter-relationships between the many dependent paciabiles of a large 
system. 

b. It provides a visual check on the character of all closed 
paths. In areal, passive system the overall effect of any feedback must 
be negative. Therefore, at each node touched by one or more loops, at 
least one of the loops must develop negative feedback. 

c. Each variable is visibly related to every other variable 
through a lattice of transmissions. Thus the effect of eliminating a 
dependent variable is readily apparent. 

d. Small sections of a large network are easily isolated and 
simulated. By application of the rules for reducing the flow graph, a 
single-path transmission can be obtained, relating the variable of in- 
terest to the forcing function, and the behavior of the isolated section 
can be predicted analytically. By comparison of the simulation results 
and the analytical solution, the validity of the simulation technique may 


be verified. 


4.3 Flow Graph of an Electrical Network. 
4.3.1 Elementary Flow Graph Theory and Operations. 


Sections 4.3.2 and 4.3.4 are largely tutorial regarding element- 
ary concepts of linear flow graph theory as applied to passive electrical 
networks. All flow graph concepts presented herein are derived from refer- 


ence [13] . For the reader already familiar with signal flow graphs, these 
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sections will serve to illustrate how the flow graph is applied in this 


study. Section 4.3.3 lists the polarity and subscripting conventions used. 


4.3.2 Relationship of the Flow Graph to the Electrical Network. 


No effort is made herein to develop the general theory of linear 
flow graphs, nor the circuit topology upon which the application of flow 
graphs to electrical circuits is based. These topics are comprehensively 
presented by S. J. Mason and Henry J. Zimmerman [13] . Only elementary 
concepts of flow graph theory are applied to the circuit of Fig. 4-1; 
these concepts are explained as they are applied. 

The signal flow graph, as it is used herein, is a graphical representa- 
tion of the set of simultaneous, linear differential equations that de- 
scribe the inter-relationships of the circuit variables. The variables 
are the nodes of the graph; they are inter-connected by transmissions 
which express the dependence of a variable upon one or more of the other 
variables. At this point it is necessary to recognize the clear dis- 
tinction between the nodes and branches of an electric circuit, and the 
nodes and transmissions of the flow graph describing the circuit. Cir- 
cuit nodes are points at which two or more circuit elements (impedances, 
admittances, generators, etc.) have a common terminal; the circuit ele- 
ments themselves comprise the branches. The nodes of the flow graph, on 
the other hand, are the circuit variables, usually current and voltage, 
which are inter-related by the graph transmissions. The functional de- 


pendence of node (variable) N, upon the remaining variables is expressed 


by 


Nj => N; We (4-1) 
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where Ni is the source node of transmission Thy. 
Flow graphs may be used equally well with both loop and nodal analy- 
sis. The latter has been chosen as best suited to the requirements of 


the model, specifically, items b. and c. of Section 4.2.1 above. 


4.3.3. Elementary Flow Graph. 


For purposes of illustration the flow graph of a segment of the cir- 
cuit of Fig. 4-1 will be developed. Fig. 4-2 shows the circuit to be 
considered; the node and circuit element designations correspond to 


those of Fig. 4-l. 


L 13 Re mM 





Figure 4-2 
Illustrative Circuit 


The polarity and subscripting convention followed through is as 
follows: 
a. Positive currents enter nodes. 
b. All node voltages are referred to the common electrical 
midpoint of the system, here shown as a common return. 
c. The voltage subscript and the first letter of current 


subscripts indicate the node at which they are measured; the second 
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Subscript on currents indicates the direction from which the current 
enters the node (a- above, b- below, 1 - left, r-right). 

In formulating the flow graph for a given system it is usually not 
necessary to express the inter-relationships between the variables alge- 
braically before diagraming them in the flow graph; indeed, it is general- 
ly redundant to do so. Such a procedure is followed in this development, 
however, so that the algebraic character of the flow graph is clearly 
established. 


The transformed equations of the circuit of Fig. 4-2 follow: 


Typ = 7 CLj-in + Ler? (4-2) 
Beer mtn 5 Couauae nour s (4-3) 
Ter = CE ~ Ex * Upp Orbe 75h, (4-4) 
Ey ca nee (Gi + R, ) in ar Len (0) a (4-5) 
Ios a le 7 Lop (4-6) 
yn = Em = Epa Wy, OEny7 Gilat ie ) ie 
=e = a Es Ce + CFOS (4-8) 
Lint = Lor a ew A Rs; (4-9) 


Note that the equations are formulated so that the functional relation- 
ship of each dependent variable is expressed explicitly once and only 
once. The flow graph of equations (4-2) through (4-9) is shown in Fig. 
4-3. | 


4.3.4 Operations on the Signal Flow Graph 
a 


The flow graph of Fig. 4-3 is a complete analytical expression, re- 
lating each of the system's dependent variables to every other and to the 
forcing functions <- the initial conditions and the independent variable, 


oa . « For computer simulation, however, further operations on the 


m 
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flow graph are desirable. 

In equation (4-5) and the corresponding transmission of the flow 
graph the term a oar Re) I yy appears; this implies a time differentiation 
of ip, (2) - Differentiators should be avoided where possible. In elec- 
tronic analog computers, internal noise causes wide voltage excursions 
across amplifiers used as differentiators. Numerical differentiation has 
inherent limitations which render its use inadvisable. hia] By absorp- 
tion of the node ar] , the transmissions from vn to er » and 
from oy to Ter can be expressed without the differentiators. 
Absorbing a signal flow graph node corresponds exactly to the elimination 
of a variable from the corresponding set of simultaneous equations. The 
governing concept in any transmission of the flow graph is, of course, 
that the total transmission between any two nodes, not absorbed in the 
transformation, must remain unaltered, 

Figure 4-4 illustrates the transformation. In Fig. 4-4(a), Es is 


absorbed. Direct transmissions from ivy are given by the ratios, 


_(s Liat R,) a uc ier R,) 
S Le (5 List FR) 


ed to each of the two. nodes adjacent to Ey which are partially dependent 


The initial condition, 1); (o)L , is shift- 
db (2 


on that variable. In Fig. 4-4(b), each of the ratios formed in the pre- 
vious operation, itself forms two transmissions - one is a gain, the 
other corresponds to a time integration. 

The direct transmissions between the nodes of the pair Cnn ) 
form a closed gain loop, abe ; similarly, a closed gain loop, — th 
exists at the node pair ( ye ee ). For reasons given in Section 5, i 
such gain loops generate instabilities in some methods of simulation. 
These loops are removed by the successive transformations shown in Fig. 


4-5. The essential character of these loops is preserved if they are 
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Figure 4-4 


Removal of Differentiators 
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Figure 4<5 


Removal of Gain Loops 
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(a) 


(b) 


(c) 
| SS K+ AX+ AxX+---- +X... 


: w(t+@+abe--- Owe 


= wel 


b= Q: 


Figure h-6 


Removal of Self-Loop 





separated into a node-to-node transmission and a self-loop. This can be 
accomplished at either of the nodes touched by the Loop. May was 
chosen as the node at which to establish both self-loops as is shown 

in Fig. 4-5(a). The self-loops at Ley are parallel transmissions 
and can therefore be added together to form the single self-loop illus- 
trated in Fig. 4-5(b). The self-loop can be replaced by a node-to-node 
transmission of the proper form. Consider the flow graph of Fig. 4-6(a). 


By splitting node “% into,’ and %" as in Fig. 4-6(b), the transmission 


Tyy is readily evaluated: 


oe Wea (4-12) 
= WwW +ax' (4-12a) 
= W/(1-2) (4-13) 
Ys y = x (4-14) 
= W (i-4) (4-14a) 
Ty = Y/w = I7(1-4) (4-15) 


This is shown in Fig. 4-6(c). 
A similar transformation of the self-loop of Fig. 4-5(b) yields the 
graph, Fig. 4-5(c), which can be used as a basic for simulating the cir- 


cuit of Fig. 4-2 on either an electronic analog or digital computer. 


4.3.5 Resume of Transformations. 

By elementary flow graph transformations, certain transmissions have 
been eliminated which would have caused difficulty in computer simulation. 
In carrying out these transformations, however, certain nodes represent- 
ing measurable system voltages or currents were necessarily eliminated. 
Similarly, new nodes, dimensionally equivalent to voltage or current, but 
not directly identifiable with a particular circuit voltage or current, 
were caused to appear. By continued repetition of the transformation 
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Processes described in the preceding sections, further reduction of the 
flow graph is possible. The extent to which such transformations are 
carried out depends upon the degree of resemblance that is desired bee 
tween the flow graph nodes and transmission and the parameters of the 
system being described. For purposes of this study the only transformae 
tions desired are those which remove differentiators or pure gain loops. 
The result is a flow graph which retains many of the circuit's node volt- 


ages and branch currents. 


4.4 System Flow Grap 





Observing the polarity and subscript conventions stated in 4.3.3 
the circuit of Fig. 4-1 is described by the flow graph, Fig. 4°7. Corres- 
pondence of the circuit with segments (a) through (3) of the flow graph 


is as follows: 


Fig. 4°7(a) SPS = starboard 
(b) Trim Pump 
(c) SPS = port 
(d) Emergency Propulsion Feeder 
(e) Drain Pump 


In order to keep the system flow graph relatively uncluttered, 
initial conditions are not shown in Fig. 4e/7. They cannot be neglected, 
however, and should be applied as in the illustrative circuit described in 
4.3. 

The circled nodes of Fig. 4-7 are those common to two or more segments 
of the flow graph. 

By processes described above, differentiators, gain loops, and selfe- 
loops are removed from the basic flow graph to yield the graph, Fig. 4-8. 
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All noneunity transmissions are labeled by one or two lower-case letters, 
and each node not otherwise identified in terms of circuit currents or 
voltages is assigned upper-case letters for reference. The transmissions 
are tabulated in Table 4-1. In Table 4-2 are listed all nodes of Fig. 4-8 
to which initial conditions are applied, together with the initial condi- 


tions. 
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Table 4-l 


Transmissions of Final Flow Graph 
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Table 4-2 


Initial conditions 


Flow graph nodes and corresponding circuit variables at t= 0, 


Node Circuit Variables 
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Circuit Values 
(Refer to Fig. 4-1) 


Component Inductance Capacitance Resistance 
Number (microhenrys) (microfarads) (ohms) 
i oe 700D8 
2 8.10521 205955 O 
3 eee 37 93 09464 594 
l 4.1518 ,032114 .000 796* 
> OO .005358 .000 796* 
6 1.01399 004 414 ~O175 
i 6.4795 0477 Ol 
8 6.9658 -OS2i2 
9. 6.7083 ~OL894 
LO 3.8406 mOL ie? 
ihe ~ 61,076 002411 
i 49.4 henrys 5.65 farads 
slag enh. ~OLOL 
14 JUNC 
Is 140. 


3 
This value is for the two diodes in their conducting state; in the 
blocking state, they would be 4,260 ohms. 
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2. System Simulation. 
5.1 Analog and Digital Computers as Simulation Devices. 
5.1.1 General. 


One of the requirements imposed upon the mathematical model in the 
previous section was that it be amenable to system simulation by computer, 
either analog or digital. This requirement has been met with the signal 
flow graph. Each transmission of the flow graph, Fig. 4-8, is of a form 
readily programmed on an electronic analog computer; similarly, each node 
is either a summing junction or signal distribution point, which may also 
be readily programmed. With regard to digital computer simulation, routines 
are available for the solution of a set of simultaneous linear differential 
equations; such equations are readily formulated in the proper format from 
the flow graph. The choice of vehicle for the simulation then is based up- 
on the relative advantages/disadvantages of the computers available with 


respect to the system to be simulated. 


5.1.2 Electronic Analog Computer. 


The advantages of simulation by analog computer are as follows: 

a. The effect of varying system parameters, initial conditions, 
and independent forcing functions is immediately observable at the output 
terminals of the operational amplifiers. 

b. The system can be readily AMeanaed or reduced with little 
effect on the overall program. 

c. The continuous and instantaneous character of analog computer 
transmissions obviates the necessity of removing gain loops from the flow 
graph. 

On the other hand certain disadvantages which occur with analog simue- 
lation are discussed in the following paragraphs. 
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The useful linear range of operational amplifiers is limited. The 
lower limit is a function of the internal noise of the amplifier, and for 
those variables being recorded, it is also a function of the sensitivity of 
the recording instrument; at the upper limit, saturation occurs. For all 
but a very simple system, these limitations generally necessitate scaling 
the variables such that their voltage analogs always remain within the 
linear range of the amplifier. Additionally, it is often desirable to 
scale the independent variable, herein considered to be time, in order to 
extend the observation interval for short duration effects, or to compress 
the interval of long term responses. The degree of difficulty of proper 
scaling increases in proportion to the increasing order of the system to be 
simulated. 

In removing the several differentiating transmission that occurred 
in the basic flow graph, a number of the system's impedance parameters lost 
their separate identity as coefficients. Each of these parameters was 
caused to appear in two or more new coefficients in different combinations 
with other parameters. On analog computers, the significant figures of the 
coefficients are normally set as fractions on potentiometers; their order 
of magnitude is programed using the gain, or attenuation, associated with 
the input and feedback impedances of the amplifiers. The accuracy to which 
the potentiometers may be set is generally limited to four significant 
figures, and is usually less. Unless each parameter is consistently repre- 
sented in each of the several coefficients where it appears, the character- 
istic roots of the system may be shifted sufficiently to appreciably alter 
the system's response. In fact, for a system which has roots with small, 
negative, real parts, small inconsistencies in the several coefficients may 


be sufficient to drive the real part of these roots positive, resulting in 
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an apparently unstable system, truly an anomaly when the system is a 
passive electrical network. This difficulty can only be minimized by 
retaining a sufficient number of figures of each parameter-factor of the 
POeErielenes to insure against error in the fourth significant figure of 
the coefficients themselves. 

The passive parameters of the circuit model established in this study 
are considered constant throughout the duration of the transient; however, 
more complete data concerning resistance vs. load characteristics of the 
siliconecontrolled rectifier, and knowledge of the exact configuration of 
the power supply circuitry may lead to a model of this device in which the 
resistance is both time and load variant. Considerable expertise may be 
required to accurately program these non-linearities on the analog computer. 

The number of operational amplifiers required for large system simula- 
tion may be prohibitive. Consider the solution by repeated integration of 


the n-th order differential equation shown below in Laplace transform no- 


tation: 
n-l ‘ 
C s a; s! 
= _ iso (5-1) 
Ss Dyess 
=o 
Now, let 
\a/ = ame IN (5-2) 
4) ‘ 
_ es 
d_ bjs 
J=o 
Then, 


N=! 


Wy i (5-3) 
}z0 


Cc 


Equation (5-3) is illustrated by the flow graph, Fig. 5-l. 
Thus the solution can be obtained by n successive integrations. In 
addition to the n amplifiers required for the integrations, a sign inversion 
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must be effected after every alternate integration to provide the proper 


signs of the feed-back and the feed-forward transmissions. These inver- 


sions may be accomplished by a single amplifier for each direction. The 
minimum number of amplifiers required is therefore N+2 , in addition to 
those that may be required to establish the proper magnitude of coefficients. 
The circuit shown in Fig. 4-l is a thirty-second order system; based upon 
the preceeding discussion at least thirty-four amplifiers would be required 
to provide a solution to the differential equation relating an input-output 
pair. Considering the flow graph, Fig. 4-8, as a basis for simulation, the 
number of amplifiers required is expanded considerably. Each of the closed 
loops that form the main lattice of the flow graph provide negative feed- 
back to all nodes touching the loop. There are, however, two integrators 
per loop and to provide negative feedback on an analog computer requires 

an odd number of amplifiers between an input quantity and the feedback. 
Thus each loop requires an inverting amplifier in addition to those required 
for integration. By placing the inverters in branches common to two loops, 
the number of these inverting amplifiers required is minimized. With such 
an arrangement, certain of the voltages and currents will be of opposite 
polarity from that shown on the flow graph, but this is of small concern. 
There are, in addition to the loops making up the lattice, nineteen other 
transmissions which will require a sign inversion in an analog simulation. 
The number of amplifiers required for system simulation, preserving as 
accessible nodes all circuit node voltages and branch currents shown in 
Fig. 4-8, is then 3/2 the number of integrators plus 19, i.e., 67. As 
with the solution of the differential equation, this number does not in- 
clude amplifiers required solely to establish the proper magnitude of co- 
efficients. The maximum number of amplifiers available at the time of 
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this study was 60, and this lack of sufficient amplifiers precluded using 
the analog computer for a full system simulation. Where sufficient ampli- 
fiers are available, or when lower ordered systems are studied, its other 
limitations notwithstanding, an analog simulation would be of considere- 


able value in a study of this type. 


5.1.3 Digital Computer. 


The discussion which follows is based upon the use of Fortran language 
with a large, general purpose computer, specifically the Control Data Corp- 
oration 1604 Computer. 

The advantages of a digital computer simulation are: 

a. System changes are readily programmed from the modified signal 
flow graph. 

b. No scaling is required. Quantities varying over many orders 
of magnitude are readily accommodated. 

c. The different time frames for the various time duration 
effects may be stipulated in the same program, subject to limitations dee 
scribed below. 

d. The impedance parameters - inductance, capacitance, and resis- 

tance =- need be specified only once; they will be retained in memory withe- 
in the computer with an accuracy of eight significant figures. 
Thereafter, all computations involving these parameters will be executed 
using consistent values of the entering arguments, and the accuracy with 
which the computer coefficients will be used far exceeds that attainable 
with the analog computer. 

e. Programs for the simulation of systems of high order can be 
handled within the large memory of the computer. 

f. System non-linearities can be programmed to whatever degree 


of accuracy may be desired. 
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ge Recorded results may include tabulated as well as graphical 
data. 

There are corresponding limitations of the digital computer which are 
discussed in the following paragraphs. 

The solution of the system's differential equations is accomplished 
by numerical integration. The largest time increment consistent with an 
accurate solution is generally considered to be 1/100 of the period of 
the fastest oscillation. If sustained high-frequency oscillations are 
super-imposed upon lower frequency oscillations, or upon long time-con- 
stant damping, many hours of computer time may be required to generate a 
steady-state solution. If, on the other hand, the damping ratio of the 
high frequency waves is relatively large, they will be attenuated after a 
few cycles and a larger integration time increment can then be used. 
Unless the approximate locations of the system's roots are known, 
several computer runs may be required to determine the sequence of time 
steps and the proper time to introduce each in order to generate the steady 
state solution in minimum computer time. 

Rhee eat ive nature of the digital computer solution forms the basis 
of the problem of pure gain loops which was mentioned in Section 4. Con- 


sider the transmission from w to z shown in Fig. 5-2. 


= wl 
$ 
| "aa cata | 
Ww v4 4 % 
Figure 5-2 


Pure Gain Loop 
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Algebraically, 


(5-4) 


ye 
it 


w- ay -fydt 
y x (5-5) 
a J 


A simultaneous solution of these three equations, for \w/(t)=!u(t) yields 


(5-6) 


the following: 


| 
3 (t) = e7 Hat (5-7) 


Now, consider the solution of the system by numerical integration for At 
very small, and as>>! . Table 5-1 shows the values of the various 
variables at time zero, and after each of the einer three time increments. 
Where used, € represents a value very small compared with the other 


quantities with which it is summed. 


Table 5-1 





Thus, when @ is large, the highest power term dominates, and the output 
increases in magnitude at each iteration by a factor approximately equal 

to Q  , and it alternates in sign. Ignoring the feedback effect of the 
integration path, it can be seen that the values of the output represent, 
in part, the series expansion 
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When Q<]} , the series can be terminated after 1 terms with an error 
less than the magnitude of the (+1) th term [13] - For Q >], however, 
the series must always include the remainder term, the reth term of equa- 
tion (5-8). The error introduced by omitting the remainder is seen to 
increase without bound as the series is expanded. When a rectangular 
integration scheme is employed, the series will expand by one term for 
each time increment; for more accurate numerical integration processes, 
such as the fourth order Runge-Kutta method flo, 18] ,» the system For- 
tran equations may be executed by the computer more than once for each 

At , causing an even more rapid expansion of the series. The two 
techniques which, in general, provide more accurate numerical solutions to 
differential equations - smaller increments of the independent variable, 
and refined estimates of the slope between steps - will, for the case of 
the gain loops, cause a more rapid divergence of the error. The problem 
can be eliminated by absorbing one of the nodes and transforming the loop 
into a node-to-node transmission as discussed in Section 4.3.4 and illus- 
trated in Fig. 4-6. 

A further drawback of digital computers as simulation devices is the 
fact that the effects of parameter variation are not readily apparent. The 
extent to which the effects may be observed depends upon factors not direct- 
ly associated with the problem to be solved. Two such factors are the access- 
ibility of the computer to the engineer, and the type of on-line display 
equipment available with the computer. 
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5.1.4 Other Simulation Schemes. 

Investigation of other simulation schemes was limited by the time 
available for the study. One simulation method worthy of consideration 
for any future study is the use of a digital analog, (hybrid) computing 
scheme, by which it may be possible to select the advantages of both 


computers and to minimize their disadvantages. 


5.2 Digital Computer Program. 
The flow graph, Fig. 4-8, shows the system to be described by an 


inter-connected network of transfer functions. Everywhere it appears, 
the Laplace transform variable, s, is in the denominator and is of first 


degree. Thus, 





Kee = K 
Kin s+P 


By cross-multiplying and rearranging terms we obtain 
sXout = K Xin — FP Xout 
and inverting this to the time domain yields the first order differential 


equation: 


a fest) = K Xin ae Xout 


Thirty two such equations were formulated fromthe system flow graph. 
These and the equations describing the inter-connections were used with a 
computer library routine, Subroutine INTEG 1 [18] to provide a numerical 


solution to the thirty-second order system. 


5.3 Simulation Problems. 

The first attempts at simulation were made using only a small section 
of the distribution system. It was with this small section that the pro- 
blem of simulating gain loops, described in Section 5.1.3, was encountered. 
Even in the small network chosen, the cause of the rapid divergence was not 
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readily apparent. The network chosen for the trial simulation was the 
circuit model of the drain pump. This network forms a part of the illustra- 
tive circuit, Fig. 4-2; its flow graph, after the removal of the differen- 
tiating transmission, is shown in Fig. 5-3. The forcing function is a 


current into the circuit node L from some other part of the system. The 





Figure 5-23 
Signal Flow Graph of Drain Pump 
closed gain loop, — be touches nodes Top and Ler . Gain trans- 
13 


missions connect these nodes with nodes § and [,,, respectively. Thus 
the divergence caused by the gain loop has almost instantaneous effect not 
only upon the nodes which touch it but also upon any node having an input 
signal which is a gain-only function of the nodes touching this loop. In 
the test circuit then, four of the six nodes were rapidly divergent. Re- 
cognition of the problem's cause occurred only when the values of all the 
variables were printed’ after each execution of the equations by the comp-~- 


uter,. 


5.3.2 Transmission Gains. 

The problem of accuracy in potentiometer settings for gains and other 
coefficients, when such values were computed values involving more than a 
single parameter, was described in Section 5.1.2. When "slide rule accuracy" 


values were used as input data for the digital computer, an instability 
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occurred there also. When the values of the parameters themselves were 
used as input data and all coefficients and transmissions computed withe 
in the program, the results indicated a stable system which agreed with 


the analytical solution. 


5.3.3. Time Increments. 

As stated in Section 3., the lumped circuit equivalent model of the 
transmission line was developed for a best frequency match of 0.1 mc. It 
was expected that this would be close to the dominant frequency of the 
system. NeGeee for all simulations of parts of the system which ine 
cluded one ar more of the transmission line T-sections, the model was 
found to be unstable for integration step sizes greater than 0.01 microe 
seconds. Many times the instability caused by too large a time step was 
not apparent until the program had progressed through a large number of 
iterations. Low frequency oscillations, of initially large magnitude, 
would be damped but then the high frequency perturbations would become 
relevant and eventually cause exponential growth. Without a complete 
determination of the system's characteristic equations, it is not poss- 


ible to more precisely predict the dominant frequencies. 
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6. Results and Conclusions. 
6.1 Results. 

A method for developing a system model and for simulating that model 
on a digital computer was developed and is shown. The simulation program 
was not perfected, and never ran correctly for the entire system. The 
time increment needed to keep the numerical integrations from diverging 
was very small; the computer program was running two hours in order to 
simulate fifty three microseconds. 

The method of development is believed to be good, and future studies 
should benefit from it. The incompatibility of gain loops and numerical 
integration methods is a fact that should be known before digital computer 


simulation is started. 


6.2 Conclusions. 


6.2.1. Routines for the solution of systems of differential equations, 





such as INTEG l, are useful tools for simulation. 


6.2.2. A valid way to test the simulation of the system or its components 
is by allowing it to change from some known initial state to another, 
known, final state; to either start with zero current and voltage and to 
"close the switch on the battery" or start with energy stored in the system 
and allow it to be dissipated. The flexibility of the digital computer 
allows a single branch or the entire network to be tested in this manner 
with only slight changes in the program. 

This check-out technique was time consuming, however. A preliminary, 
much shorter method of indicating convergence, was to put the expected 
final values in the program as initial conditions and see if they were 


stable. 
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6.2.3 Although it was plausible to call the transmission lines lossless, 
including a slight resistance in the model for them would have been a good 
idea. If nothing else, computation time would be decreased because of 


the smaller transient decay times. 
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APPENDIX B 
DISTRIBUTED PARAMETERS of CABLING 
B.1l Inductance. 
B.1.1 Special Considerations. 

Although numerous reference works on transmission line theory are 
available, two geometrical factors of the submarine d-c distribution sys- 
tem require special attention in the development of the distributed in- 
ductance relationships: (1) the location of the multi-cable system close 
to the high permeability steel hull, resulting in fields which vary azimuth- 
ally as well as radially; (2) the proximity effect - when the distance 
of separation of adjacent cables is not large relative to the linear cross- 
section dimension, the adjacent sides of the cables will be linked by more 


flux than will the far sides. 


B.1.2. Method of Images. 


Consider the magnetic field established by a long, straight, cylin- 
drical, current-carrying conductor located in a region of low magnetic 
permeability and far removed from other conductors and from any region of 
higher permeability. The constant flux lines of the magnetic field in the 
region are a series of circles concentric with the axis of the conductor. 
If the conductor is now brought to the vicinity of other conductors, and 
near a boundary separating the low permeability (air) region from one of 
high permeability (steel), the flux lines are no longer symmetric in 
azimuth about the center of the conductor; the total flux distribution 
is influenced by the relative polarities and the magnitudes of all the 
currents, and by the shape and proximity of the air-steel boundary. 


For the case where all the currents are parallel to a plane, vertical 
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boundary, ! the boundary conditions of the horizontal and tangential 
components of the magnetic flux lines on the air side of the boundary 
are identically satisfied as follows:~ 

a. For each actual conductor postulate an image conductor 
located on the opposite side of the boundary, and at equal distance from 
the boundary as the actual conductor, i.e., at the image point of the 
original conductor. 

b. Establish in each image conductor a current equal to that 


in the corresponding actual conductor, multiplied by a factor of 


Ky = cere (B-1) 


where Mic and MU, are respectively the permeabilities of the steel and 
air. 

c. Consider the region surrounding the conductors, both actual 
and image, to have permeability Ma. 

Thus, the problem of determining the strength of a single, highly 
asymmetric field reduces to the relatively straightforward superposition 
of 2N symmetric fields, where " is the number of actual conductors. 

When the permeability of the second medium is very much greater than that 
of the first, as is the case with steel and air, Ky approaches unity, and 
both actual and image currents may then be considered equal. 

1 

The submarine hull is actually not a plane surface at any point, but 
rather is circular in corss-section. However, its center of curvature is 


sufficiently far from the centers of the cables that it may be regarded as 
a plane for the purposes of this discussion. 


No attempt is made herein to justify the method of images as it is 
applied to magnetic fields. For the theory of magnetic images see [1] 
and [8] . 
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A typical configuration of submarine d-c power cables and steel hull 
boundary is shown in Fig. Bela. The equivalent system, from an inductance 
point of view, with image conductors added and steel hull removed, is shown 
in Fig. B-lb. All subsequent calculations for the distributed inductance 
are based upon the model of Fig. B-lb, which may be modified to include 
more, or fewer cables, and different relative locations of the opposite 


polarities. 


B.1.3. Flux Linkages. The magnetic flux linkage of each conductor is the 


sum of; 
a. Self-linkage of a current by the flux established by that 
current 
(1) Partial flux linkages by lines of flux within the 
conductor; 


(2) Total linkage by all lines outside the conductor; 
b. Mutual-linkage of each current with flux established by the 

current of each other conductor (See Fig. B-2). 

(1) Partial linkage of current in conductor 1 by flux of 
current in conductor 2 at digtance pP from the center of conductor 2, 
where 

D-k< 9 < Dth 

(2) Total linkage of all current in conductor 1 by the flux 
of conductor 2 when fP = D+ te ; 

The flux linkages of a.(1) and b.(1) above will be modified by the 
"skin effect" of high frequencies. In the self-linkage term, this effect 
will be taken into account by adjusting the computed dec value according to 
tabulated ratios. For the mutual flux linkage term, the skin effect is 
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Figure Bela 


() Typical d-c cable configuration 


for hull mounted power cables, 


Figure B-lb 
Image model of Fig. B-la. Hull has been removed and image conductors 


added. Note: the image conductor has the same polarity as the corres- 
ponding actual conductor. 


rs 





integrally treated in the general development of the appropriate formula. 


Bo1.4. Internal Self-linkage. The basic d-c relationship and the skin 


effect receive considerable attention in transmission line reference 
works; consequently only the necessary results are presented herein. 


For the d-c case 


ALT: 
8 


For skin effect calculations, a table of internal self-inductance ratios 


d ii(internal) = webers/meter (B-2) 


is available [12] ; the entering argument for the table is the quantity 
(i = ele = (B-3) 

where [2 is the radius of the conductor, 

uJ) is the radian frequency, 

A is the permeability of the conductor, 
and P is the resistivity of the conductor. 
The table may be used for values of x & 100. Where x falls between two 
values for which the inductance ratios are listed, linear interpolation 
between the two given values will yield results to a sufficient degree 
of accuracy for most calculations. It is of interest to note the wide 
variation of the ratio: 
at x = 0, L,./Ly.= 1.000; at x = 100, Ly. /Ly. = 0.028. For values of 


x > 100, a close approximation to the inductance ratio is given by 


bas % 
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B.1.5. External Self-linkage. The desired effect of introducing image 
currents was to obtain azimuthal symmetry of the flux field of each con- 
ductor considered separately. The contribution to external linkage of 

each conductor by its own current, when such symmetry is established, is 


given by the expression 


A ii(external) = 2 x 107/ 


ln 1/r, webers/meter (B-5) 
It is important to note at this point that the above relationship can be 
derived only when all conducting paths, forward and return, of the trans- 
mission circuit are considered. Any development of inductance or flux 
linkage relationships, for flux external to a conductor, will result in an 
indeterminate mathematical expression if the complete circuit is not 
considered. For example, the analytical expression for the external 
self-inductance of a single conductor, without reference to any other 
conductor, is 

°° 


t-> 
L, (external) =2x 107/ ln r | henrys/meter (B-6) 
rT = 


r 
Cc 


B.1.6. Partial Mutual Linkages. The necessity of considering the proxi- 


mity effect is generally obviated by assuming that the center-to-center 
separation of the conductors is much greater than the radius of the con- 
ductor. In the submarine dec distribution system, such an assumption is 
unwarranted; the separation of the cables is generally of the order of 
only 2 - 3 times the conductor radius. A general development is pursued 
wherein the conductor area of the cable is an annular ring of outside 
radius r. ‘and inside radius ry . (Fig. B-2). Radius r, is a function 
of frequency and at dec becomes zero; at high frequencies it is approxi- 


mated by the difference between the conductor radius r,. and the skin 
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penetration depth § , As frequency is decreased however, the error of 
this approximation becomes very large. For the same entering argument 

x that was used with the internal selfeinductance skin effect corrections, 
tabulations of aec/d-c resistance ratios for solid, round wires are avail- 
able. Now, if the resistance is considered to be constant across the annu- 
lar conduction area, and the resistivity of the conductor invariant with 


frequency, then 


FRRac x 1 Ge = rs) = 
ine = The (B-7) 


Solving equation (B-7) for rs 


r = r [t= Rye / Rae meters (B-8) 


For x > 100, the skin depth ¢, for copper is given by 


-2 
arc = 6.62 X10 meters (B-9) 


TF 


and the inside radius of the annular ring is therefore 


r = pe & meters (B-10) 


The magnetic field intensity of conductor 2 at radius eP is given by 


B = Mole | webers/meter* (B-11) 


= 5 P 
and encloses that portion of the annular ring of conductor 1 which is 


shown shaded. To find the area of this shaded segment proceed as follows: 


a, Find the area A, of the triangle 0, 0, P, 
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A, = = D E San CS meters ¢ (B-12) 


b. Find the Area As of sector 0, Pe S 


A. = ¢ ae. meters (B-13) 


c. Find the area A3 of sector 0, iv re 


A; = + ts Xe meters” (B-14) 


d. Obtain the area A, bounded by the arcs Pe S and P, T,, 


and by the straight line st. 


Ag = Ag aes (A, is Ay ) 


A, = & Lrt«e -(Dr sinx,-P'A)] 15) 


meters© 
where 
z 2 ae 2 
% = cos "(2 5a") — @-16) 
<i c 
2 2 pe 
GB. = cos "(+ +#_—_ ay E ) (B-17) 


e. In similar fashion it will be found that the ares Ag ene 


closed by the arcs Pe S and P, Ty and the straight Line ST. is given 


by 
As m= > [ nS x, aarti SH ae eae. )] (B-18) 
xX; = COS Se ) @-19) 
3 Sh; Oe 2 aes 


80 





f. The area Ag of interest is twice the difference between 


the areas Ay and Ag ‘ 


Bye = <2 (Ay a As ) (B-21) 


With the area determined, the current linked by the flux of conductor 


2 can be determined, and the flux linkage computed in the usual manner: 
= po As 


where 


iL 
D> 
1 |® 
if 
ES 


Cha = UE) 2 FP 


D+ 
= Wi, 1 ; 
ro) <i 272 (r2- 2 rz fA 2 dp (B-22) 


(Note A, is included within the integral because of its dependence on). 


B.1l./. External mutual linkages. For the same reasons outlined in Sec- 
tion 5 above, the development of external flux linkage relationships can- 
not be pursued without considering all current carrying paths of the cir- 
cuit. However, the result may be expressed as a summation of terms of the 
form 


N : _ hy | (B-23) 





| 
Mm 
> 
S 
O 
a 
> 





B.1.8. Conclusions and Results. 

The relationships derived in the preceeding section were programed 
for digital computer calculation. Input data included cable spacing, 
number of cables of each polarity, their arrangement, and cable diameter. 
In addition, Table 15-2, reference 12 was read in. 

Table B-2 shows the variation of inductance with frequency for one 
cable type, the largest cable in the distribution system. The maximum 
error from neglecting the proximity effect is shown to be approximatee- 
ly 3.5% at 10 kc. Since the effect increases with increasing conductor 
Size, this is the worst case. 

The values of inductance in Table B-l were computed using the proxi- 
mity effect, but in view of the above results, neglegible error would be 


involved if this effect were ignored. 
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Cable Type: 


No. of Pairs: 





Frequency 


0 


100 kc 


1 Mc 


10 Mc 


100 kc 


1 Mc 


10 Mc 


La. /L 100 kc 


TABLE B-1 


EQUIVALENT DISTRIBUTED 


INDUCTANCES AND RESISTANCES 


SSGA 800 SSGA 800 SSGA 650 
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SSGA 400 


1 


- 19204 
~ 17392 
. 16800 


«15885 


- 12360 
2.82058 
8.81258 


27.8678 


3 2 3 
Distributed Inductance, 
© 39246 - 27196 «41250 
- 34301 « 23793 - 36339 
31875 22220 «33535 
© 29277 « 20550 ~ 31285 

Distributed Resistance, 

~02057 - 03086 0.2589 
» 66219 - 99328 » 14139 
2.07430 3.11378 2.32869 
6.55855 9.84664 7.3640 
1.14 1.14 1.13 


1.10 


SSGA 300 


1 


- 20192 
- 18339 
~ 17990 


- 16876 


017255 
3.33856 
10.41232 


32.9266 


1.10 


DSGA 125 


1 


«15358 
- 14009 
«13583 
«12999 


« 205407 
2.5928 
8.0331 


25.4029 


1.10 





TABLE B-2 


The Variation of Distributed Inductance with Fregvency 





a. Considering skin and proximity effects 
b. Considering skin effect alone 


(2 pairs of SSGA 800 calbes, separated 1.6 inches) 


Frequency Inductance, a Inductance, b 
cps Microhenrys/meter 
0 0.27196115 0.27187966 
10 , 27219310 27173755 
100 275593862 . 26287288 
10° . 25767800 . 22602490 
10° 24721412 . 21185031 
10° . 23793139 . 20727483 
10° 22225427 . 20587376 
107 . 20550048 20542194 
108 . 20535731 , 20527907 
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B.2. Distributed Capacitance. 


B.2.1. Single Conductor Cable. 


Enclosure (5) to reference [2] provides a tabulation of shunt capa- 
citances for SSGA type cables, sizes 3 through 2000. For the multi- 
conductor case, where half the cables carry current in one direction and 
the other half carry the return current, the total distributed capaci- 
tance is the tabulated value multiplied by the number of conductor pairs 
(parallel addition of the capacitances of like polarity cables) and divided 
by two (the series addition of capacitance between cables of opposite polar- 
ities). An additional factor is included to convert from microfarads-per- 


thousand-feet to microfarads-per-meter. 


2.2 Two Conductor Cable. 
The relationship for the capacitance between two conductors enclosed 
within a conducting sheat is developed by Atwood [1] . After adjusting 


for units, the equation becomes: 
a Inf2s. e=s4 AU 
a rt +57 
ks relative permittivity of the dielectric surrounding 


the conductors 


S = displacement of the center of the conductors from the 
center of the cable 


Y% = radius of the conductor 
Y) = radius of the enclosing sheath. 


For a dielectric with relative permittivity k = oe 


ISSR 
—— nfs Bese = 


iz re o eG 
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B.3. Distributed Resistance. 

The d-c resistance per unit length was first calculated for a single 
conductor for each of the cable sizes. This value was then adjusted for 
the various cable configurations in the same manner as the capacitances: 
a three pair cable being considered as three parallel sets of two series 
resistances. 

Using the methods of reference [12 ] , these dec values were adjusted 
to account for skin effect, and 0.1 mc values were found. The increase 
was, in general, a resistance 25 to 30 times the d-c value. The results 
are summarised in Table Bel. 

The resistance calculations were done in conjunction with the induct- 
ances discussed in section B.1, although the resistance does not include 


a proximity effect. 


B.4. Calculation of Lumped Equivalent Circuit Values. 


The preceeding sections of this appendix derive the equivalent ine 
ductance, capacitance, and resistance per unit length of the various cable 
configurations. The definitions in Section 3.2 and Equations (3-30) are 
the means of computing the equivalent circuits values from them; the line's 
length is the other factor. For cable number 3, at 0.1 me, 

Q = 67.056 meters 
L. = 0.363389 microhenrys/meter 
G = 0.00256 864 microfarads/meter 


Vo = === = 32,7312 x 10° M/sec 
LC 
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Table B.4. shows the values used in this model, the characteristic 


impedances, and the loss parameters (which are subsequently neglected). 


y_- RQ | 
ak oF, 
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APPENDIX C 
Passive Circuit Models 
C.l. The d-c Motors. 

Figure (C-1) is an equivalent circuit of a separately excited d-c 
motor of the type driving the Trim and Drain Pumps. The known data, from 
which an entirely electrical circuit is derived, include: 

the output power 

the line voltage 

the motor speed 

the resistance and the time constants of the 


field and armature paths 
the moment of inertia. 


a CIRYA 
ee 
— > 
i 
Ip 
L 
Vie : ae f 
{- 
Re Te ap — 
"6 
Figure C-1 


Equivalent Circuit of Separately Excited d-c Motor 
The dynamic behavior of the motor is illustrated by means of a linear 


signal flow graph, Fig. (C-2). 


Vi | Vin Vo Le Im | aby 





+ 


Ke Ky 
so Td Ww Em 
Figure C-2 
Flow Graph of Motor's Dynamic Behavior 
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In the ensuring development, all relations are based upon a consistent 
system of units: when nameplate data is used to obtain quantitative re- 
sults, additional conversion constants may be required. The use of con- 
stants - in lieu of linear functions of the field current, [f¢ - as propor- 


tionality factors relating armature current, I., to developed torque, !j  , 


a? 
and relating motor speed, © to generated emf, Ea is justified when the 
time constant of the motor's field is compared with the period of the expect- 
ed transients. The relatively large time constant of the field means an 


almost invariant current during the transient. 


Now consider the electrical circuit of Fig. (C-3): 


Lea Ra 
j 





Co 


Figure C-3 
Equivalent Circuit 


The two circuits, the motor and the R-L-C, are similar. Viewed from 
an external point, they are equivalent (under the assumption of constant 


field current) if: 


ee aii 

Yh Ym 

| — Ke Ky 
SCy + Vpn sotf 


| 
= (sT/K, Kv) +(f/ Ki Kw) 


Cr = Cn Ry= Ke Kv 
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7 
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The values of Ke, Av and f may be derived from the given machine data 


and from an assumed overall conversion efficiency, Ne 


” & Power delivered to load = alee 
Power supplied to input terminals Fy (C-1) 
Using (C-1) and the assumed efficiency, the power input is found: 
re P. /% (C-1la) 
Then: t= eC / Vv (c-2) 
and: Ip = V./ Re (C-3) 


(At steady state, all the terms with time derivative dependence will 
go to zero.) 
ee (c-4) 
The power dissipated as heat in field and armature resistance losses 


is 


2 
Rtec, = [pe +T Ra (C-5) 


The difference between power supplied and electrical losses is the 


mechanical power developed at the shaft: 


eons =P - Rec. (C-6) 
But 


Hnech S10) “CU \ca2) 


and an = Ki Le (C-8) 


Eliminating ae from (C-7) and (C-8) and solving for K+: 


Kes . (c-9) 
a 
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The friction coefficient may now be found using the relation 


Td =Ki-La =wf (ree 
f - Ke La /w (C-10a) 


The voltage, Eneas, generated by the armature is a function of the coil 
geometry, the rotor speed, and the magnitude of the stationary magnetic field. 


The coil geometry is fixed and the stationary magnetic field is constant. 


Therefore, 
= = Ky: CW) (C-11) 
Ky = Em /w (C-11la) 
but ore = V; - ila Ro (C-12) 
therefore ocr = (V;- nen Ra )/w (C-13) 


With the Trim and Drain Pump, the equivalent circuit values are found by 


Pi= 125/.90 = 139 = 103.8 kw 


Jy 
Ly 


(103.8 kw 250 kv)= 415a 


250/61.7=4.05a 


i 


*- Ly = 415-4.05=4lla 


Prech= 103.8 -[(4-05)*61.7) +(411)"(0.0175)] x10 





= 99.9 kw 
~ 93.9 x10? x 60 _ 0.664 newi-m 
Ke = Se Sepcioteaaror = amp 
= ‘ v 4.) 1x 6GO0~«1O = Oeran -newt-m 
- Za Xe On 5 KO aa rad. 
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RR = ©.664 x O 6644 as 0.5384 2 
L— 0.741 = 


Lp = 0.8* 61.7 = 49.4 hy 
La = 0.0175x0.035= 0.612 millihenrys 


Cp = 59x 4.214*10% ~s cep 
(.664) 


C.2. The Battery. 


The storage battery is extremely nonlinear. Its internal resistance 
varies with the current, the rate of change of current, and with the degree 
of charge fille It has inductance [15] , and two qualities which resemble 
capacitance -overvoltage and polarization. 

Figure C-4 is the model advocated by reference [15] and [19] for short 


time studies. 





Figure C-4 
Equivalent Circuit for Battery 


The model is derived in [19] to explain the short time (one second) 
effects, under relatively low current, seen when the applied voltage is switch- 
ed. There it is asserted that for a single cell, C will be about 0.016 
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farads per rated ampere hour and that Ro will be 0.2 ohms per ampere hour. 
L is shown to approach 0.15 microhenrys for a single cell of the rating 
considered in this simulation, although the effects of the inductance are 
not visible in the reference's results. 

Reference | 14 | discusses the overvoltage that appears while the battery 
is being charged and is nearly at its rated capacity. The overvoltage's 
decay is also shown, under open circuit conditions and it appears to have 
a time constant in the neighborhood of 65 seconds, with the overvoltage 
starting at 22.5% of the nominal cell voltage. This reference also gives, 
as experimental data, the value of 40 micro-ohms per cell under short cir- 
cuit conditions of less than 0.5 seconds duration. Extrapolating the values 
of Ro and C from reference 19 to the 6,500 ampere-hour battery gives an 
overvoltage decay time constant of 135 seconds, which is considered a good 
comparison to the 65 seconds from reference [214 | ° 

The battery in this study will be 126 cells, each of 2.02 volts, in 
series. | 

Since this study is concerned alee transients, one case should concern 
the presence of overvoltage. For this, the following values would be used: 

E.: 255 v (nominal battery voltage) 
va 62.4 v (initial overvoltage) 
L: 975 microhenrys 

Ca: 60 kfarads 


R,: 0.00392 ohms 
R2: 0.00108 ohms 


This network and initial voltage match the open circuit overvoltage 
decay and the steady state internal resistance. 

In the initial simulations, the overvoltage is not considered and the 
battery is considered to be discharging. For this situation, the capacitor 


is removed from the circuit and the two resistances are combined. 
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C.3. The Voltage Regulator. 


When power is flowing from the a-c side (through a high voltage dec 
stage) to the d-c side, the active elements are SCR-1 and Diode D-2 [7 | ; no 
current passes through D-l or SCR-2. Paths 1 and 2, in Fig. C-5 are alter- 
nately followed as SCR-1 is switched to its open or blocking state. Where 


single elements are shown, there are actually parallel strings of series 


———=  — ——— _ /——— s(—( 
| 3 +450, 
SCR-1\7 7\ 0-1 
a-— — —@ — — < tft 
Ve ‘ 
y ~--—+-<-@-—--; 
/ 
LOAD D-2 |7N \/ SCR-2 
; 
\ a ; ; 
a de a a 
Figure C-5 
Current paths, a-c to d-c operation 
elements. 


Reference [7] gives a value, not otherwise available, from the charact- 
eristic curve for a single, open, SCR. This value of E/I (for the combination) 
is taken as defining a resistance, and this resistance is taken as the linear 
approximation to the SCR in its conducting state. Figure (C-6)illustrates 


this approximation. 
Ze 


2s PPROXIMATION 
SLOPES 






OS 0a 
Figu re (C-6) 


SCR approximated by 
a resistance 


R= 3.03 milli-cohms 


Figure C-6 
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For the diode, only some power dissipation figures were given. For the 
first approximation, it is assumed that power is dissipated only during -the 
conduction period, and that the diode current has a square waveform. This 
leads to a value of 0.796 milliohms for the diode in its conducting state. 

When either the diode or -the SCR is in its blocking state, it has 
practically all the 450 volts across it. Without data as to the "reverse 
characteristics" of such devices, a first approximation is to consider them 
as perfect switches, with an infinite impedance in their reverse direction. 
A second approximation is to consider them as passing a leakage current 
equal to one-ten thousandth of the normal operating current. This order of 
magnitude is based on that of the small solid state devices. On this basis, 
the blocking resistance of the Diode and of the SCR is taken to be 4,260 
ohms. 

The results is that two linear circuits are alternately switched into 


the network, as shown in Fig. (C-7). 


0 450y 
3.03«10° ae AE 
eae G 
: rh Ff 
: -6 
4260 4.510 
Figure C-7 


Equivalent Circuit of dec Side of Voltage Regulator 
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